liz22_del_14_33
seascape <- seascape_probability(reefoutline=reef_map, habitat=benthic_map)
liz22_del_14_33 <- st_read("/Users/rof011/coralseed/code/2021_2022_HighPRT_sim_examples/run_day_12038_liz22_del_14_33/run_day_12038_liz22_del_14_33.json")
## Reading layer `run_day_12038_liz22_del_14_33' from data source
## `/Users/rof011/coralseed/code/2021_2022_HighPRT_sim_examples/run_day_12038_liz22_del_14_33/run_day_12038_liz22_del_14_33.json'
## using driver `GeoJSON'
## Simple feature collection with 605000 features and 3 fields
## Geometry type: POINT
## Dimension: XYZ
## Bounding box: xmin: 145.4446 ymin: -14.71132 xmax: 145.4806 ymax: -14.69039
## z_range: zmin: 3.55 zmax: 3.55
## Geodetic CRS: WGS 84
particles <- seed_particles(input=liz22_del_14_33, seascape=seascape, probability="additive",
limit_time=24, competency.function = "exponential",
simulate.mortality = "typeI", simulate.mortality.n = 0.1,
return.plot=TRUE, silent=FALSE)

## ## coralseed model summary: ####
## Importing 5000 particle tracks
## [No seed set, random draws used]
##
## Time start = 2022-12-18
## Time end = 2022-12-19 00:57:00
## Total dispersaltime (hrs) = 1.03958333333333
##
## Competency at t6 = 639 / 5000 larvae
## Competency at t12 = 1182 / 5000 larvae
## Competency at t24 = 2100 / 5000 larvae
##
## Survivorship curve
## Mortality at t6 = 19 / 5000 larvae
## Mortality at t12 = 120 / 5000 larvae
settlers <- settle_particles(particles, probability="additive", silent=FALSE, return.plot = TRUE)
##
## 2061 / 5000 larvae settled
##
map_coralseed(seed_particles=particles, settle_particles=settlers,
seascape_probability=seascape, subsample=500, restoration.plot=c(100,100),
show.tracks = TRUE)